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ABSTRACT: 

Since its original publication in 1916 under the title " The algebraic theory of modular systems^\ 
the book by F.S. Macaulay has attracted a lot of scientists with a view towards pure matcmatics 
(D. Eisenbud,...) or applications to control theory (U. Oberst,...). However, a carefuU examination 
of the quotations clearly shows that people had only a look at the first three chapters respectively 
dealing with the resultant, the resolvent and the general properties of ideals but mostly did not 
look at the last chapter dealing with the so-called "inverse system" or, to be fair, only to the 
particular example of finite dimensional vector spaces on the residue field of a local ring. 

The basic intuitive idea is the well known parallel existing between ideals in polynomial rings 
and systems of PD equations in one unknown with constant coefficients. Accordingly, it becomes 
evident that Grobner bases should be everywhere ... under the condition to understand what is 
inside (we invite the reader to have a look at the book in order to discover the difficulty of such 
a task !). As a first step done in the last paperback edition, a kind of glossary has been exhibited 
between Macaulay language/definitions and modern corresponding concepts. However, nothing 
has been done for the last (and most difficult !) chapter and the purpose of this paper will be first 
to correct this gap. Roughly speaking ( and though striking it could be !) this chapter is based on 
the (correct) descripton of " systems" in the sense given by the american school of PD equations in 
the seventies (H. Goldschmidt, D.C. Spencer, D.G. Quillen,...) and we shall prove that the heart 
of the chapter is "i/ie use of the Spencer operator on sections'" {not on solutions, despite formal 
power series are used !). Of course, this result, (not known to our knowledge) will give by itself 
quite a new insight on the chapter. 

The next idea will be to use modern "algeraic analysis" also pioneered in the seventies (M. 
Kashiwara, B. Malgrange, V.P. Palamodov,...) in order to extend these ideas to general differen- 
tial modules and systems of PD equations. Again, we shall prove that the use of Groebner bases is 
central ... but well hidden as we shall prove that one cannot avoid Spencer cohomology, involution 
and delta-regular coordinates in order to understand in an intrinsic way the local computations 
done by Macaulay through the so-called " localization^^ technique. 

Hints will finally be given to generalize these results to the variable coefficients case and the 
corresponding control identifiability problem. Many explicit examples will be provided for illus- 
trating the main results that provide new hints for applying computer algebra to algebraic analysis. 

KEY WORDS: Partial differential equations, Macaulay inverse system, algebraic nalysis, com- 
mutative algebra, homological algebra, localization, duality, computer algebra, Grobner bases. 

1) INTRODUCTION: 

With only a slight abuse of language, one can say that the birth of the formal theory of systems 
of ordinary differential (OD) equations or partial differential (PD) equations is coming from the 
work of M. Janet in 1920 [9] along algebraic ideas brought by D. Hilbert at the same time in his 
study of sygyzies for finitely generated modules over polynomial rings. The work of Janet has then 
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been used (without any quotation !) by J.F. Ritt when he created differential algebra around 1930 
[21], namely when he became able to add the word " differentiar in front of most of the classical 
concepts concerned with algebraic equations, successively passing from OD algebraic equations to 
PD algebraic equations. In 1965 [5] B. Buchberger invented Grobner bases, named in honor of his 
Phd advisor W. Grobner, whose earlier work done in 1940 on polynomial ideals and PD equations 
with constant coefficients provided a source of inspiration [8]. However, the approaches of Janet 
and Buchberger/Grobner both suffer from the same lack of intrinsicness as they highly depend on 
the ordering of the n independent variables and derivatives of the m unknowns. 

Meanwhile, commutative algebra, namely the study of modules over rings, was facing a very 
subtle problem, the resolution of which led to the modern but difficult homological algebra with 
sequences and diagrams. Roughly, one can say that the problem was essentially to study properties 
of finitely generated modules not depending on the "presentation" of these modules by means of 
generators and relations. This very hard step^^ is based on homological/cohomological methods like 
the so-called " extensiovi^ modules which cannot therefore be avoided. 

In order to sketch this problem, let us present two simple examples. In the first case with 
standard notations, everybody will understand at once that integrating the second order OD 
equation y = with m = n = 1 is equivalent to integrating the system of two first order 
OD equations — = 0, if' — 0. However, even with m = n — 2 and the same two un- 
knowns u,v in both cases, it is not evident at all that integrating the second order PD equation 
di2U — d22V — u = is equivalent to integrating the system of two fourth order PD equations 
dii22U — dx222V - d22V - u — 0, c?iii2M — ^1122" " duu = (cxercisc !). 

As before, using now rings of " differential operators^^ instead of polynomial rings led to dif- 
ferential modules and to the challenge of adding the word " differential in front of concepts of 
commutative algebra. Accordingly, not only one needs properties not depending on the presenta- 
tion as we just explained but also properties not depending on the coordinate system as it becomes 
clear from any application to mathematical or engineering physics where tensors and exterior forms 
are always to be met like in the space-time formulation of electromagnetism. Unhappily, no one of 
the previous techniques for OD or PD equations could work !. 

By chance, the intrinsic study of systems of OD or PD equations has been pioneered in a totally 
independent way by D. C. Spencer and collaborators after 1960 [23], using jet theory and diagram 
chasing in order to relate differential properties of the equations to algebraic properties of their 
symbol, a technique superseding the " leading term" approah of Janet or Grobner but quite poorly 
known by the mathematical community, even today. 

Accordingly, it was another challenge to unify the "purely differentiaF approach of Spencer 
with the "purely algebraic" approach of commutative algebra, having in mind the necessity to use 
the previous homological algebraic results in this new framework. This sophisticated mixture of 
differential geometry and homological algebra, now called " algebraic analysis" , has been achieved 
after 1970 by V. P. Palamodov for the constant coefficient case [16], then by M. Kashiwara [11] 
and B. Malgrange [14] for the variable coefficient case. 

It is only in 1990, thanks to the work of U. Oberst, that such a theory has been applied with 
success to control theory [15]. Then the things went on rather fast towards computer algebra 
and many packages now exist for computing the extension modules and related concepts (See 
for example: http: / /wwwb. math. rwth-aachen.de / OreModules') . However, of course, many difficult 
problems are left and we provide details about a few of them, having in mind the recent workshops 
on Grobner bases and applications successively held at RISC in 2006 and 2008. 

When a given system of linear PD equations of order q is given, it defines by residue a dif- 
ferential module M over the underlying ring D of differential operators. Then it becomes today 
possible to decide by means of computer algebra the class to which M belongs among n classes 
ranging from free, torsion-free, reflexive, ... , to projective and free. However, the set of elements 
of M , namely of finite linear combinations of the unknowns and their derivatives modulo the given 
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PD equations and their derivatives, such that each of them does satisfy at least one PD equation 
for itself, provides the torsion submodule t{M) and M is torsion-frec if its torsion submodule is 
zero. An open but useful problem, independent of the previous classification, is now to classify 
elements in t{M) by means of the nested chain of n differential submodules: 

= i„(M) C i„_i(M) C ... C ti{M) C toiM) = t{M) C M 

A first question is thus to determine the classes and the "gaps" in the above chain, as indeed, in 
many known explicit situations, a few intermediate modules do coincide. A second question is to 
avoid the double extension modules in the exact sequences: 

— > tr{M) — > tr-i(M) — > ext'jj{ext'jj{M,D),D) 
generalizing the well known case r = 0, namely the exact sequence: 

— > t{M) — > M -U homD{homD{M, D), D) 

where the morphism e is efined by e{m){f) = f{m),\lm € M,'^f G homD{M, D). A new technique 
proposed for such a purpose is to bring back the system of order g to a system of order one without 
any equation of order zero, called " Spencer form" , which generalizes for PD equations the Kalman 
form existing for linear OD systems in control theory. The interest is to provide new domains of 
applications and we sketch the underlying idea on a simple academic example. 

Again, with evident notations and m = l,n = q = 2, the system ^222/ = O^duy = clearly 
determines a torsion module t{M) = M. The torsion elements z' = d\y and z" = d2y do not satisfy 
the same decoupling PD equations as z' only satisfies d2z' = while z" satisfies d2z" = 0, diz" = 
and we have the nested chain with strict inclusions: 

= t2{M) c ii(M) c to{M) = t{M) = M 

the classification being obtained through the dimension d or rather codimension cd = n — d oi the 
differential modules generated by the respective torsion elements as wc have indeed cd{Dz") = 2 
and cd{Dz') = 1. More generally, if m e M and cd{Dm) — r, then m e tr-i{M) but m ^ tr{M), 
according to the Hilbert-Serre theorem stating that cd{M) = cd{V) for ay differential module 
defined by a system of PD equation with characteristic variety V. Of course, the same decoupling 
type problem can be asked for any engineering system in gasdynamics or magnetohydrodynamics 
(MHD) but we do not know a single work existing towards such a classification in view of the 
difficulty of the mathematical framework involved. As an ultimate goal, a particularly important 
problem should be to study the dependence of the previous classification on parameters when the 
system depends on certain constant parameters, a result generalizing the controllability problem 
for OD systems in control theory where n = 1 only !. 

Accordingly, the hope should be to have a computer algebra package providing the classes and 
eventually generating elements. A particularly interesting case should be to characterize r-pure 
modules, namely modules M such that there exists an integer < r < n with tr{M) = and 
tr-i{M) = M. Equivalently, M is r-pure whenever cd{Dm) = r^m S M. Surprisingly, for 
constant coefficient systems in one unknown, such a concept has been discovered in 1916 by F. S. 
Macauly in the last chapter " The inverse system and modular equations" of his famous book " The 
algebraic theory of modular systems" but this result is still not understood today [13]. 

In fact, while looking at this last hapter since many years, we were convinced that the double 
picture of p. 67 was nothing else than, sise by side, the matrix of the "system" in the american 
sense, combined with an ordering of the various possible (formal) " truncated sections" , again in 
the american sense, the underlying idea being to cancel successively the terms of order zero, then 
zero and one, ... and so on. However, it is only a few months ago that we suddenly understood the 
true reason for supposing, as a crucial assumption indeed though it is only presented as a purely 
technical argument (p. 89), that the module under study was "unmixed'. We explain thereafter 
this point by translating the " old' language into the " new" one. 
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First of all, what Macaulay was calling "module" is now called "ideal' ([13], glossary of the 
last edition). Accordingly, the properties (prime, primary, unmixed,...) attributed to an ideal o in 
the ring A = k[x] = k[xi, Xn] of polynomials in the n indeterminates xij Xn with coefficients 
in the field k are now, along with the modern setting of commutative algebra, attributed to the 
residual module A/ a. Then we got in mind that, in the study of an r-pure differential modules, a 
delicate though expected theorem is stating that the corresponding characteristic variety is " un- 
mixed' too, with the same meaning as above, that is the underlying algebraic set is the union 
of irreducible components of the same codimension r. This analogy was thus giving rise to the 
challenge of relating the work of Macaulay on unmixed polynomial ideals to the study of purity 
for differential modules. In particular, the extension to this new framework of a " localization cri- 
terion" , provided by Macaulay in the classical setting, constitutes one of the main results of this 
paper and provides new hints for applying computer algebra to algebraic analysis. Most of the 
concepts and results presented in this paper are new and illustrated by explicit examples. 

2) MOTIVATIONS 

We start presenting a few motivations from commutative algebra, then apply localization theory 
to systems of OD equations and finally generalize the results so far obtained to systems of PD 
equations. 

Let k be field of characteristic zero and x — (xii ■■■,Xn) be indeterminates over k. We introduce 
the ring A = k[xi, Xn] of polynomials with coefficients in k and various classes of ideals. The 
set of maximum ideals is denoted by max (A) with elements m,..., the set of (proper) prime ideals 
is denoted as usual by spec{A) with elements p,... and the set of primary ideals with elements q 
such that ab G q,b ^ q ^ a G p = rad{q) that is a'' G q for a certain integer r € N. The importance 
of primary ideals lies in the fact, largely emphasized by Macaulay, that any ideal a can be written 
as a finite irredundant intersection a = qi Pi ... fl q^ of primary ideals, called primary decomposition. 
Setting pi = rad{qi), we obtain at once the prime decomposition rad{a) = pi n ... n p., though 
sometimes this new decomposition may not be irredundant with strict inclusion p^ C pj for certain 
couples of indices (i, j). In this case one uses to say that the component defined by pj is embedded 
into the component defined by pi in the algebraic set defined by a. Also, for any prime ideal p, 
we denote by cd{A/p) = n — d{A/p) the codimension of A/p with d{A/p) = trd{Q{A/p)/k) the 
transcendence degree of the algebraic extension of the field of fractions of the integral domain A/p 
over the field k. For an arbitrary ideal a, the codimension is usually denoting the minimum among 
the codimensions of the components defined by the minimum prime ideals in the corresponding 
prime decomposition, which are therefore not embedded. 

DEFINITION 2.1: An ideal a C A is unmixed if cd{A/ pi) = ... = cd{A/ps) in a primary decom- 
position and we have therefore pi ^ pj,V(i,j). Otherwise a is said to be mixed. 

We now present a few examples that will be used in the sequel with a totally different approach. 

EXAMPLE 2.2: q = ((xa)', XiXs - X2) ^ rad{q) = p = (xa, X2) ^ cd{A/q) = 2. Similarly a = 
(X1.X2X3) = (xi,X2) n (xi,X3) is unmixed with cd{A/a) = 2 but a = ((Xl)^XlX2,XlX3,X2X3) = 
(Xi)X2) n (XI1X3) 1^ (xi7X27X3)^ is mixed with two minimum prime ideals and one embedded 
component. More generally, any ideal having a basis containing as many polynomials as the codi- 
mension of the corresponding residual module has been called "ideal of the principal class" by 
Macaulay who proved that any such ideal is unmixed ([13], 48, p 40,49). For a modern approach 
through regular sequences, see ([12] , VI, 3, p 183). 

EXAMPLE 2.3: = ((X2)^XlX2) = (X2) n (xi,X2)^ = qi n qz is mixed with q^ = pi c p2 = 
rad{q2) = (xi,X2) G maa;(Q[xi, X2])- 

EXERCISE 2.4: ([13], §42, p 44) Prove that a = ((xl)^ (X2)^ ((Xi)' + (X2)')X4 + X1X2X3) is 
mixed with s = 4 and only one minimum prime (Hint: use the fact that (xiX2)^(Xi) X2, X3) X4) C a 
but X1X2 ^ a)- 

The main idea is then to " transfer" the properties of an ideal C A to the residue module M = A/a 
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over A. 



DEFINITION 2.5: A module M is said to be prime {primary) if ax = 0,0 ^ x € M ^ aM = 
{a^M = for a certain integer r) though people sometimes add the prefix "co". 

Now, having in mind the so-called Chinese remainder theorem ([12], p 41), any primary decom- 
position gives rise to a monomorphism — > Qi © •■• © Qs with primary modules Qi = A/(\i and 
epimorphisms M ^ Qi — > 0, Vi = 1, s. Conversely, looking for such a situation allows to exhibit 
a primary decomposition for reducible modules (see [18], p 110 for more details). 

It is now tempting, and this too was a key idea of Macaulay, to introduce n commuting deriva- 
tives di, ...jdn for which k should be a field of constants and to introduce the ring D = k[d] = 
k[di, dn] of differential operators with coefficients in k. As D and A are isomorphic by di Xi) 
any (nonlinear) ideal of A gives rise to a system of OD/PD (linear) equations in one unknown 
only and conversely. It thus remains to use techniques for PD equations in order to study ideals or 
modules. However, the situation for a differential field K with subfield of constants k and/or sys- 
tems of PD equations for many unknowns escapes from the previous approach and we conjectured 
that they could be treated " by their owri^ , the specific situation considered by Macaulay being just 
a particular case of the general theory that we shall present in this paper. 

First of all, we sketch the technique of " localization!^ in the case of OD equations, comparing 
to the situation met in classical control theory where n = 1 and the dimension can therefore only 
be or 1. For this, setting as usual d = di = d/dt = dot, we may introduce (formal) unknowns 
y^, and set Dy = Dy^ + ... + Dy"^ ~ D™. If we have a given system $ = of OD equations of 
order q, a basic question in control theory is to decide whether the control system is " controllable" 
or not. It is not our purpose to discuss here about such a question (see [18,19] for more details) 
but we just want to state the final formal test in terms of a property of the differential module 
M = Dy/D^. Care must be taken that in the sequel, for simplicity and unless needed, we shall 
not always put a " bar" on the residual image of y in the canonical projection Dy — > M. We explain 
our goal on an example. 

EXAMPLE 2.6: With to = 3 and a constant parameter a, we consider the first order system 
= y^ — ay^ — y^ = 0, = y^ — y'^ + y^ = 0. Any engineer should want to apply Laplace 
transform y{s) — e''*y{t)dt to this system. However, using the integration by part formula 
f'^e^*y{t)dt = [e''*?/(t)]g° — sy{s) we should eventually need to know y{0) though the Kalman test 
of controllability is purely formal as it only deals with ranks of matrices [10]. Since a long time 
we had in mind that setting y(0) = was not the right way and that Laplace transform could be 
superseded by another purely formal technique. For this, let us replace "formally" d by the purely 
algebraic symbol x whenever it appears and obtain the system of linear equations : 



but we could have adopted a different choice for the only arbitrary unknown. At this step there 
are only two possibilities : 

•a 7^ 0, 1 =>no " simplification" may occur and, getting rid of the common denominator, we get an 
algebraic parametrization leading to a differential parametrization as follows: 

= Xix + a)z, = xix + l)^, = (X^ - a)z ^ y^ = z + az,y'^ = z + z,y^ = z - az 

•a = Oora=l=^>a "simplification" may occur. For example, with a = 0, setting z = y^ — we 
get x-z = that is to say i = 0. 

Recapitulating, we discover that a control system is controllable iff one cannot get any au- 
tonomous element satisfying an OD equation by itself. For understanding such a result in an 
algebraic manner, let M be a module over an integral domain A containing 1. A subset 5 C A is 
called a multiplicative subset if 1 G 5 and Vs, t € S => st G S. Moreover, we shall need/use the Ore 



Xy^ - ay^ - xy^ = 0, ly^ - xy^ + X2/^ = 



Xix + a) 3 
o y 
X -a 




x(x + i) 3 
X -a 
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condition on S and A, namely aS r\ sA ^ ^,ya £ A, s € S. 

DEFINITION 2.7: S-'^A = {s-'^a\s € S,a € A/ r^} with s-^a ~ t-'^b <^ 3u,w e A,us = vt G 
S, ua = vb. 

Next, for any module M over A, we define S^'^M = S'~^A(g)^M and ts{M) = {x & M \ 3s e 
S, sx = 0} in the exact sequence ts{M) ^ M — » S~^M where the last morphism is x — » l~^x. 

EXAMPLE 2.8: S = A - {0} ^ Q{A) = K field of fractions of A and we introduce 

the torsion submodule ts{M) = to{M) = t{M) = {x & M \ ^ a & A,ax = Q} oi M . Also, if 
p G spec{A) and S = A — ^, one uses to set S~^M = rup. 

PROPOSITION 2.9: When M is finitely generated and t{M) = 0, from the inclusion M C 
K^aM, we deduce that there exists a finitely generated free module F with M c F. 

REMARK 2.10: Though the above proposition provides a parametrization for any n in the case 
of a torsion-free module, in the particular case n = 1 there is an isomorphism M ~ t{M)(BM/t{M) 
not so well known in OD control theory. Indeed the projection onto the second factor is the canon- 
ical projection onto the torsion- free module M/t{M) which is a free and thus projective module 
when n = 1, a result allowing to split the short exact sequence — > t{M) M ^ M/t{M) — »■ 0. 
This is not evident at all on Example 1.6 and even on the very simple example — = 0. 

The comparison with Example 1.6 needs no comment at least when n = 1 and controllability 
must have to do with t{M) = when n > 2 though it is only quite later on in the paper that 
we shall be able to generalize the result expressed by the above remark. Also the extension of 
the above results to the non-commutative case D = K[d] where -ft" is a differential field with n 
commuting derivations di, can be achieved but is much more delicate ([18,19]). 

EXERCISE 2.11: When a = a{t) in Example 1.6, we let the reader prove that the controllability 
condition is now the Ricatti inequality d-|-a^ — a^Oina coherent way with the constant coefficient 
case already considered. 

Taking into account the works of Janet and Spencer, the study of systems of PD equations 
cannot be achieved without understanding involution and we now explain this concept. For 
this, let — {^i, iin) be a multi-index with length \ ^ \— + ... + We set n + li = 
(/xi, + 1,/Zj+i, and we say that n is of class i if /xi = ... = Hi-i = 0,/Kj ^ 0. 

Accordingly, any operator P = a^dfj, € D acts on the unknowns y'' for A; = 1, m as we may set 
d/^y'' = with j/q = y'' and introduce the jet coordinates yq = {y^ \ k = l,...,m;0 <| fi \< q}. 
It follows that, if a system of PD equations can be written in the form = aJJ'^y^ = with 
a £ K,we may introduce the filtred differential module M = Dy/D^ but we notice that the work 
of Macaulay only covers the case m= 1. Then we define the (formal) prolongation of with re- 
spect to di to be di^'^ = aY'y'^+i^ +^iC'k^y^ ^^'^ induce maps di : M ^ M : y^ ^ v'pi+u by residue. 

Changing linearly the derivations if necessary, we may successively solve the maximum number 
of equations with respect to the jets of class n, class (n — 1)...., class 1. At each order, a certain 
number of jets called principal {pri) can therefore be expressed by means of the other jets called 
parametric {par). Moreover, for each equation of order q and class i, di, ...,di are called multi- 
plicative while di+i,...,rf„ are called nonmultiplicative and are nonmultiplicative for all 
the remaining equations of order < q — 1. 

DEFINITION 2.12: The system is said to be involutive if each prolongation with respect to 
a nonmultiplicative derivation is a linear combination of prolongations with respect to the multi- 
plicative ones. Using Spencer cohomology, one can prove that such a definition is in fact intrinsic 
[7,17,18,23]. 

EXAMPLE 2.13: ([13], §38, p 40 where one can find the first intuition of formal integrability !) 
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The ideal q — ((xi)^, X1X3 — X2) provides the system j/n = 0, ^13 — 2/2 = which is not involutive. 
Effecting the permutation (1,2,3) (3,2, 1), we get the new system 1/33 = 0, yi3 — 2/2 = 0. As 
diysa - ^3(2/13 - 2/2) = 2/23 and diy23 - ^2(2/13 - 2/2) = 2/22, the new system 2/33 = 0, 2/23 = 0, 2/22 = 
0, 2/i3~2/2 = is involutive with 1 equation of class 3, 2 equations of class 2 and 1 equation of class 1. 

APPLICATION 2.14: t{M) = M iff the number of equations of class n is m. Otherwise there 
is a strict inclusion t{M) C M. 

PROPOSITION 2.15: ([17,23]) The following recipe will allow to bring an involutive system 
of order q to an equivalent (isomorphic modules) involutive system of order 1 with no zero order 
equations called Spencer form: 

1) Use all parametric jets up to order q as new unknowns. 

2) Make one prolongation. 

3) Substitute the new unknowns. 

PROPOSITION 2.16: For such a system, defining the character a\ = m— number of equations 
of class i, we have a\ > al > ... > a" > and only the first nonzero character is intrinsic, coming 
from the Hilbert polynomial of the module/system. 

REMARK 2.17: cdoiM) > r ^ a^"''+^ = ... = a'l = 0. In that case, we shall say that we have 
full class n,..., full class (n-r+1). Thus r = 2 in the above example. 

3) PURITY 

As a first basic fact, quite important for the study of the noncommutative case, one must carefully 
distinguish between an ideal/system and its symbol part, namely the top order part of order q 
when the system is involutive. The following example will explain the difficulty involved, hidden 
in the use of Grobner bases which are not intrinsically defined. 

EXAMPLE 3.1: The primary ideal corresponding to the involutive system of the preceding ex- 
ample of Macaulay has radical p = (x3!X2)- However, the symbol part is the homogeneous ideal 
= ((X3)^,X2X3, (X2)^XlX3) = {(X2)^X3) n (xi,X2,X3)^ with the same radical but it is a pure 
coincidence. Nevertheless, a and q have the same (co) dimension according to the famous Hilbert- 
Serre theorem that wc shall recall in Proposition 3.6. The importance of involution has not been 
pointed out clearly in the study of Grobner bases. However, it is clear that the two polynomials 
(Xs)^ and XiXa ^ X2 generate q in the previous example but the corresponding homogeneous ideal 
at order 2 is now ((X3)^)XiX3) with radical (xs) providing a wrong dimension. 

DEFINITION 3.2: The characteristic variety V of an involutive system of order q is the alge- 
braic set defined by the radical (care !) of the polynomial ideal in K[x] generated by the m x m 
minors of the characteristic matrix {a]J^XiJi) where \ iJb\= q. Of course, when m = 1 we recover the 
radical of the ideal we started with. 

If m e M, then the differential submodule Dm C M is defined by a system of OD/PD equa- 
tions for one unknown only and we may look for its codimension cdD{Dm). In the commutative 
case, looking at the annihilators, wc get ariri{M) C ann{Dm). In particular, if M is primary its 
annihilator is a primary ideal q with radical p and we have q C ann{Dm) C p,Vm e M as a 
possible characterisation. Accordingly, if M is prime, then ann{Dm) — p,Vm G M. 

EXAMPLE 3.3: In Example 2.13, with the primary ideal q, then 2/2 and 2/3 are killed by p though 
y is killed by q. The situation changes completely with the corresponding homogeneous ideal as 
2/1 is killed by (x3, (X2)^) and 2/12 is killed by (X3,X2) though 2/3 is killed by (xi,X2,X3)- 

Even in the noncommutative case of systems with coefficients in a differential field K, one can 
prove with the homological techniques of algebraic analysis ([2,18]) : 

PROPOSITION 3.4: tr{M) = {m & M \ cd{Dm) > r} is the greatest differential submodule of 
M having codimension > r. 
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PROPOSITION 3.5: tr{M) does not depend on the presentation of the module as it can be 

defined inductively by the exact sequences : 

— > tr{M) — > U-iiM) — > exe^{ext'jj{M,D),D) 
if we start from t-i{M) = M and to{M) = t{M) when r = 0. 

Thanks to (http://wwwb.math.rwth-aachen.de/OreModules), this proposition is essential for 
the use of computer algebra and allows to refer to the Spencer form. In fact, the situation is exactly 
similar to that of control theory with the Kalman form. 

PROPOSITION 3.6: cdoiM) = cd{V) = r O a'i'"" ^ 0, a5'"''+^ = ... = = <^ ^r-i(M) = 
... = to(Af) = t(M) = M. 

We may therefore define as in the Introduction: 

DEFINITION 3.7: M is v-pure O tr-i{M) = M,tr{M) = ^ cd{Dm) ^ r,\/m & M. In 
particular, M is 0-pure iff t{M) = 0. 

The following key results using a kind of partial localization^^ generalize the similar ones first 
obtained by Macaulay ([13], §82) and provide a technical test linking purity and involution. From 
now on we shall only consider the constant coefhcient situation, considering xi, •■■,X«-r just like 
parameters ([13], §77, p 86), but most of the results can be extended to the variable coefficient 
situation, though with a lot of work more. 

MAIN THEOREM 3.8: If cdoiM) = r one has the exact sequence: 

tr{M) M ^ k{xi, X«-r) ® M 

Proof. According to the definition of involution, the system made by the PD equations of class 
l+...+class (n—r) is also involutive for di, dn-r and thus also for di, d„ by adopting the order- 
ing {dn-r+i, dmdi, dn-r)- It allows to define a differential module and an epimorphism 
Mr — > Af ^ as M is defined by more equations. Now, as cdoiM) = r, we have tr-i{M) — M and 
each torsion element of t{M) — M surely satisfies at least r PD equations involving successively 
dn, dn-r+i- As for the other equations, they should only include di, dn-r and this is just the 
way to construct t{Mr) by considering the exact sequence: 

> t{Mr) > Mr k{xi, ...,Xn-r) ® Mr 

As it is well known that localization preserves exactness [3,22], this exact sequence projects onto 
the desired one. 

Q.E.D. 

REMARK 3.9: The above technique allows to generalize for arbitrary m the condition obtained 
by Macaulay for m = 1 ([13], §41, p 43 and §43, p 45) that we translate into modern language as 
another proof (See [3], IV, §1, exercise 9). In fact, if a C A = k[xi, X«] is such that cd{A/a) = r, 
then a is unmixed (A/a is r-pure) if and only if ^(xi, Xn-r)P G a =^ P £ a. For if o = qifl.-.nqs 
with, say cd{A/qi) > r, then 3S € qi and P ^ qi, P G q2 n ... n qs so that SP S a does not require 
Pea. Conversely, if no such qi exists and SP G a, then P G a : {S) = o. Indeed and more 
generally, if b G A is an ideal and a C a : b 7^ a, then b C pi = rad{qi) for some i. To prove this, if 
a : b — c and b is not contained in any pi, then one can find 6 G b, 6 ^ pi, Vi and c G c, c ^ o with 
6c G a C qi and thus c G qi,Vi that is c G a and a contradiction. 

COROLLARY 3.10: ([13], §41, p 43) M is r-pure ^ Af -> k{xi, ...,Xr) ® M is exact. 

COROLLARY 3.11: M is r-pure 30 ^ M ^ L with projective dimension pd]j{L) = r if 
r > 1. 
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Proof. (Compare to [2], p 494 and [18], p 553) As M is r-pure, then tr{M) = that is Mr is torsion 
free ad we may use k{xii Xn-r)(E)Mr in order to parametrize Mr exactly as we did for embedding 
a torsion-free module into a free module. According to Proposition 2.16, the parametrization now 
depends on a""'' arbitrary unknowns z, that is we may embed Mr into a""'' copies of k[di, dn-r] 
when coming back to the differential framework. After substitution into the original equations, 
the equations of class 1 up to class (n — r) disappears for the z as they are automatically satified 
by the parametrization. The number of nonmuliplicative derivatives is < r — 1 (care) for each of 
the remaining equations of class (n — r + 1) up to class n. But such a number is just the way to 
know about the projective/free dimension by constructing a resolution of M or a Janet sequence 
for the system ([17], p 146). 

REMARK 3.12: In actual practice it is important to notice that the partial localization "kills'^ 
the PD equations of class 1 up to class (n — r — 1) {care again !) because of the compat- 
ibility conditions provided by the involutive assumption. Moreover we now understand why 
Macaulay ([13], §79, p 89) was always dealing with unmixed ideals o or pure modules A/a. It 
is known ([2] and [18], Proposition 3.173, p 549) that an A-module M is r-pure if and only if 
cd{A/p) = cd{M) = r for any p G spec{A) appearing in the prime decomposition of rad{annA{Mj) 
and no embedded primary components occur in a primary decomposition of annA{M) or equiv- 
alently ass{M) = {p g spec{A) \ 3x € M, ann{x) = p}. Example 2.3 shows that the second 
condition is needed. Accordingly, any prime or primary module is pure. 

EXAMPLE 3.13: A; = Q, to = 1, n = 3, g = 2, r = 1. The module M defined by the involutive 
second order system = 2/33 = 0, (f>^ = 1/23 = 0, (f>^ = 2/13 = is not pure. Among the compati- 
bility conditions we have d2^^ — di^^ = As only the class 3 is fuU, the localization is done by 
tensoring with k{xiTX2) and we get $^ = ^'^^ ■ Also, as yi and 2/2 arc killed by ^3 they are in 
<o(M) = M but not in ti{M). Hence there is a "gap" because ti{M) = t2{M) and t/3 generates 
t2{M) as it is killed by (di, ^2,^3) and t^{M) = by definition. In order to get a first order pre- 
sentation (though not a Spencer form) we may introduce z^ = y, z'^ = yi, z^ ~ y2, z^ = y^ and Mi 
is defined by 3 equations of class 2 and 2 equations of class 1 (exercise). With respect to (^1,^2), 
Ml is not torsion- free but not a torsion module as t{M{) is generated by z'^ as already noticed. 
Finally we have ass{M) = {(xs), (xi) X2, X3)}as another way to check that M is not pure. 

EXAMPLE 3.14: fc = Q, m = 3, n = 4, g = 1, r = 1. The module M defined by the first order in- 
volutive system y\ = 0, 2/| = 0, y\ = 0, = is 1-pure. We notice that Mi is defined by the 
only divergence- free condition and is thus torsion- free. Indeed, tensoring by A;(xi,X2;X3) iii order 
to localize, we get the parametrization = —^y^ — = —X2Z^ — ilP' = Xs-^^^i?/^ = Xs-^^ 
and we have an embedding M (Z L with L generated by (z^, z^) satisfying only z\ = 0, z| = 0. Ac- 
cordingly, L admits a resolution O^I?'^^£)^^L— >0 with morphism (Pi, P2) — * (-Pi '^4, Pid^) 
and pd{L) = 1. 

4) INVERSE SYSTEMS 

Let X be a differential field with subfield of constants k = cst{K). The ring D = K[d] is filtred 

by the order q of operators and we have K = Dq C Di C ... C -Doo = D. Accordingly, the module 
M is filtred by the order q of the linear combinations yq = Dgy allowing to describe elements of 
M and we have the inductive limit Mq C Mi C ... C M, C ... C Mqo = M with djMg C M^+i. 

DEFINITION 4.1: We define the system R = homxiM, K) = M* and set Rq = homxiMq, K) = 
M* as the system of order q in order to have now the projective limit R = R^o ii, ^ 

... Ri Rq. Accordingly, if a system of PD equations of order q is given as before, then 
fq&Rq'.y^^f^GK with a^'^f^ = defines a section at order q and we set foo = f & R for a 
formal solution. 

REMARK 4.2: In the case of an involutive system of order q in solved form, the matrix (aJJ^) and 
the corresponding prolongations for increasing ordersin allow to express the principal jets from the 
parametric jets (called " complete set of remainders" by Macaulay) and exacty describe the upper 
part of the picture drawn by Macaulay in ([13], §59, p 67 and §68, p 79). Similarly, for any q >0 
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the following commutative and exact diagram: 
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allows to define the symbol of Rq and the upper row is again exactly describing the lower part of 
the same picture through the use of truncated formal power series. The symbols are a modern way 
to describe the " compartments" of Macaulay. The use of a basis (1, 0, ....), (0, 1, ....) and so on for 
the parametric jets, after ordering them, brings a diagonal matrix as explained by Macaulay ([13], 
§59, p 67). 

DEFINITION 4.3: A ''modular equation" E = a^f^ = of order q is just a "practical' way to 
write down a section /g G i?g if <| |< g by using " effective sections" and "formal coefficients" 
instead of " effective coefficients" and "formal unknowns" in the case of PD equations. 



The following proposition generalizes the results of Macaulay to arbitrary systems with variable 
coefficients ([2], [18]). 

PROPOSITION 4.4: When M is a £>-module, then ii is a £»-module too. 

Proof It is clear that D, as an algebra, is generated hy K = Dq and T = Di/Dq with Di = K(BT. 

Let us define: 



{af){m) = af{m) = f{am) Va £ K,ym € M 

(^/)(m) = ^/(m) - /(^m) = a'di eTym^M 

It is easy to check that dia = adi + dia and that $,ri ~ rjS, = [^, rj] is the standard bracket of vector 
fields. We finally get {difY^ — {d'if){y^) — difj^ — fj^^i. that is exactly the Spencer operator. In 
particular didj = djdi,\li,j = l,...,n (exercise) and diRq+i C Rg leading to diR C R though 
decreasing the order, contrary to what happens for M. 

Alternatively and in a coherent way with differential geometry, if we have a linear system ^'^ = 
defining Rq and its first prolongation $^ = 0,di^'^ ~ defining Rq+i, as already exhibited in 
Section 2, a section /g+i e Rq+i over fq e i?, satisfies both al^f^ = and a^'^/^+i. +dial'^fl^ = 
as equalities in K with <\ n \< q. Applying di to the first and substracting the second, we get 
O'l'^i^ifll - f^+u) = 0- Accordingly, we obtain: 

fq+, G Rq+i ^ difq+, £ Rq ^ E = alf^ = ^ d^E ^ a^^{dj^ - = 0,V/ G i? 

but diE is of order q with <| /x |< g whenever E is of order q + 1 with <\ 1/ \< q + 1. When 
K = k we recognize exactly (up to sign) the operator of Macaulay ([13], §60, p 69). 

Q.E.D. 



EXAMPLE 4.5: fc = Q, m = 1, n = 1, g = 2. For the system - y ^ 0, if we set: 

/' = (1, 0, 1, 0, ...) ^ E' = a° + fill + ... = and /" = (0, 1, 0, 1, ...) E' = + + ... = 

and use {/', /"} as a basis of R over k, we have dif" = f or cquivalently diE" = E'. 

REMARK 4.6: If cd{M) = r and K = k, then q^"''+^ = 0, = and a partial local- 

ization brings the system to a finite type (zero symbol at high order) system in dn-r+i, dn) 
over the field fc(xi, Xn-r)- Accordingly, there is a finite number of linearly independent cor- 
responding sections of the localized system and thus an equal finite number of (dual) modular 
equations as in the previous example ([13], §79, p 88). In the situation K = k, we have also 
annniM) = annoiR)- Indeed, as a representative of any element of M can be written as a 



10 



finite linear combination of parametric jets with coefficients in fc, we have M C M** and thus 
ann{M) C ann{M*) C ann{M**) C ann{M) => ann{M) = ann{M*). This result generalizes the 
one of Macaulay ([13], §61, p 70) obtained when m = 1. Indeed, a^^ffi = => a'^^f^+v = by 
prolongation and thus E = a^f^ = => a'^'^d^E = a'^'^a'^fn+v = 0. 

As we already know, any primary decomposition, say with two components for simplicity, gives 
rise to a monomorphism ^ M ^ Q' (B Q" where Q' , Q" are primary modules, both with two 
epimorphisms M — > Q' ^ 0, M ^ Q" 0, respectively induced by the localization morphisms 
M Mpi,M Mp« when M is pure (unmixed annihilator) with ass{M) = Setting 
R' = homK{Q',K),R" = homK{Q" ,K) and using the fact that D* is injective, we get an epi- 
morphism R' ® R" — > — > both with two monomorphisms ^ — > i?, — > R" — > R proving 
that R' , R", R' + R", R' n R" are subsystems of R. The following proposition, not evident at first 
sight, explains the aim of Macaulay ([13], end of §79, p 89) and allows to use various subsystems 
for studying R instead of decomposing M. 

PROPOSITION 4.7: R = R' + R". 

Proof We have the well known short exact sequence ^ R' (1 R" ^ R' ® R" ^ R' -\- R" ^ Q 
where the last morphism is (/', /") — > (/' — /"). Composing the cpimorphism with the monomor- 
phism ^ R' + R" R and using the fact that the composite morphism R' ® R" ^ R is an 
epimorphism, it follows that the previous monomorphism is also an epimorphism and thus an iso- 
morphism, though in general R' fl R" ^ 0, unless we have p' + p" = A, a. situation always met with 
max{A) c spec{A) in the case of modules over a ring A which is also a finitely generated algebra 
over a field k. 

Q.E.D. 

We finally recall in a self-contained way a few results on the so-called socle and top of a module 
M over a commutative noetherian integral domain A with unit 1 ([1],[6]). First of all, we quote 
the following theorem on associated primes where both isolated and embedded components are 
needed ([3], IV, §1, exercise 11). 

THEOREM 4.8: If M is a finitely generated A-module, the sequence -» M ®peass(M)-^p 
is exact. 

Proof. If the sequence is not exact, let N be the kernel of the morphism on the right. If ass{M) = 
{pi,...,ps}, let us consider the defining exact sequences ^ Ni M ^ Mp.,\/i = l,...,s. By 
definition, we have N = (iNi and it is well known that N ^ ass{N) ^ 0. In that case, let 
p e ass{N) c ass{M), that is to say p = pi for a certain 1 < i < s. Again by definition, one can 
find X € N c Ni such that p = pj = ann{x). But x € Ni ^ 3s j £ Si = A — pi, SiX = because 
of localization and thus a contradiction. One could also say that ann{x) C pi for some i whenever 
7^ X € A'' c M as it is well known that Upi is the set of zero divisors for M. But x G Ni and we 
conclude as above. 

Q.E.D. 

REMARK 4.9: Using this theorem while taking into account Remark 2.10 and Proposition 4.7, 
we are able to relate controllability to the irreducibility of a module ([12], p 179 and [18], p 110) or 
equivalcntly to the study of maximal subsystems corresponding to the torsion-free specialization 
in control theory. Such a result will be the heart of this section when dealing with the top of a 
system and seems totally new. 

EXAMPLE 4.10: Let A = Q[x,y,z],M = A/a with a = {x'^,xy,xz,yz) = pi n p2 n where 
pi = {x,y),p2 = {x,z) are the two minimal primes (isolated components of the characteristic 
variey) and m = {x,y,z) e m,ax{A) (embedded component). Then ass{M) = {pi,p2,tn} where pi 
kills z, p2 kills y and m kills x. It follows that x belongs to the kernel of M — *• Mp^ © Mp^ and 
cannot be killed by any s G A — m. 

Keeping in mind the bricks needed in order to construct a house, a basic idea in module theory is 
to look for the greatest semi-simple submodule of a given module. For this, if m G max{A)riass{M), 
then one can find a finite number of elements x,y,... £ M killed by m. Accordingly, the map 
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X : A ^ M : a —* ax has kernel m and A/m ~ Ax C Af is a simple module, like Ay which 
may eventually be different and so on. The direct sum Ax (B Ay (B ■■■ is called the socle of M at 
m and denoted by socm{M). The various simple components are called isotypical as they are all 
isomorphic to A/m. 

DEFINITION 4.11: The socle of M is soc{M) = ®soCm{M) for m G max{A) n ass{M). It is 
the largest semi-simple submodule of M. 

We notice that the double condition on the direct sum is essential as we need not only a sub- 
module (m G ass{M)) but also a simple module (nx e max{A) C spec{A)). Also, if S', S" are two 
simple submodules of M, then M/S' © S" is the fiber sum of M/S' and M/S" over M ([12], p 88) 
but the resulting construction is not natural and will provide a motivation for duality in order to 
use Proposition 4.7. Finally, M is semi-simple if M = soc{M) and soc{M) = if M has no simple 
submodule, like the Z-module Z. In the previous Example S = Ax ~ A/m is the only simple 
submodule of M. 

EXAMPLE 4.12: If A = Q[x,y] and M = A/a with a = {x^,y'^,xy), then both y and are 
killed by m = (x, y). It follows that soc{M) = soCm{M) = Ay(BAx^ has two isotypical components 
isomorphic to A/m. 

LEMMA 4.13: Any niorphism f : M ^ N induces a morphism / : soc{M) — > soc{N). In 
particular, if M' is a submodule of M, then soc{M') = M' D soc{M). 

Proof. The lemma follows at once from the Schur lemma saying that, when f ^ 0, then M simple 
/ injective, N simple / surjective. 

Q.E.D. 

DEFINITION 4.14: The radical of a module M is the submodule rad{M) which is the inter- 
section of all the maximum proper submodules of M. If rad{M) = 0, for example if M is simple, 
we say that M has no radical. If M has no proper maximum submodule, then rad{M) = M. 

LEMMA 4.15: rad{M) is the intersection of all the kernels of the nonzero morphisms M S 
where 5 is a simple module. 

Proof. From the Schur lemma, the above morphism is an opimorphism and we may introduce 
the defining short exact sequences 0^-/V— >M— ^S— »0. Let us consider the exact sequence 
DN ^ M — > (BS. The image of the morphism on the right is a submodule of a semi-simple 
module and thus a semi-simple module too, which is even a direct summand. Accordingly, restrict- 
ing the choice of the simple modules in order to have an irredundant intersection, the morphism 
on the right thus becomes an epimorphism providing the short exact sequence of the lemma. 

Q.E.D. 

DEFINITION 4.16: The top of the module M is the semi-simple module defined by the short 
exact sequence rad{M) M top{M) —> 0. It can also be defined as the largest quotient 
of M that is a direct sum of simple modules. 

We have the following three useful lemmas: 

LEMMA 4.17: Any morphism f : M —>■ N induces a morphism / = rad(M) rad(N). 

Proof. Let S be a simple module. For any morphism g : N ^ S, the composition g o f : M ^ S 

vanishes on rad{M) and thus g vanishes on f{rad{M)), that is f{rad{M)) C rad{N). 

Q.E.D. 

LEMMA 4.18: If M 7^ is finitely generated, then rad{M) ^ M. 

Proof From noetherian arguments, M always contains a maximum proper submodule. 

Q.E.D. 

LEMMA 4.19: (Nakayama) Let M be a finitely generated module and A'' a submodule of rad{M). 
If L C M is such that L + N = M, then L = M. 
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Proof: Let us suppose that L ^ M. Then, from noetherian arguments again, L is contained in a 
maximum proper submodule L'. It foUows that N + L C rad{M) + L' ^ L' and a contradiction. 

Q.E.D. 

We are now ready to provide the achievement of this paper while explaining ([14], §77,79,82). 
5) MACAULAY ' S SECRET 

The crucial idea of Macaulay has been to replace soc{M) by top{R) by means of duality theory, 
in order to use Nakayama's lemma for finding generating sections (formal solutions) of R. We 
proceed in a few successive steps for working with differential modules in an effective way and 
treating the following specific examples. 

5.1) The first basic procedure is to check that M is r-pure. For this we must determine r by 
exhibiting an involutive system. As already noticed, a linear change of derivations may be needed 
in order to chock involution. The partial localization will then be used in order to check the purity 
and to deal only with maximal ideals because a prime ideal is maximum if and only if its residue 
integral domain is zero dimensional. 

EXAMPLE 5.1.1: If a = (xi,X2X3) = (Xi,X2)n(xi,X3) the corresponding system ?yi = 0,2/ii = 
0, yi2 = 0, yi3 = 0, 2/23 = is not involutive and the change xi -* X3, X2 ^ X2, X3 X2 - Xi Pro- 
vides the involutive system in solved form 7/33 = 0, 2/23 = 0, 1/22 — yi2 = 0, yis = 0, 7/3 = 0. We have 
M C Q(xi) €5 M and the localized module has the two associated maximum ideals rrii = (^3,(^2) 
and m2 = (^3,^2 - Xi) with tui + m2 = Q{xi)[d2,d3] as d2 - {d2 - Xi) = Xi- 

EXAMPLE 5.1.2: a = ((xi)^) (X2)^, X1X2) is primary because rad{a) = (xi:X2) = m. With 
D = Q[di,d — 2], the homogeneous second order system R3 = {2/222 = 0, 1/122 = 0, 2/112 = 0, 2/111 = 
0, 2/22 = 0, 2/12 = 0} is trivially involutive because its symbol is zero. The corresponding module is 
primary and 2-pure. No localization is needed and dimq{R) = dimq{R3) = 4. 

EXAMPLE 5.1.3: a = ((X3)^ X2X3 - (Xi)^, (X2)^) is primary because rad(a) = (xi,X2,X3) = m. 
With now D = Q[di,d2,d3], the homogeneous system R2 — {2/33 — 0, 2/23 — yii = 0, 2/22 = 0} is not 
involutive (see [18, p 321 for another similar example) but its prolongation R4 is trivially involutive 
with zero symbol and dimq{g3) = 1. No localization is needed and dimq{R) = dimq{R2) = 8. 

5.2) The idea is now to adapt to modules an argument already used in Remark 3.9 for ideals. If 
a = nqj is a primary decomposition with pi = rad{qi), then a : b 7^ a => bj C pj for a certain i. In 
particular, a : m ^ a for m € max{A) =^ m = p, for a certain i. 

Let us set A"^ = F as a free module and consider the short exact sequence 0— >J^F^M— >0 
where / is the so-called module of equations and let a e ^ be an ideal. We want to prove the fol- 
lowing lemma: 

LEMMA 5.2.1: / C/ :a = J7^/^aCpjfora certain i. 

Proof. Let us consider a primary decomposition I = fl/j in F and pass to the quotient by introduc- 
ing short exact sequences ^ li F ^ Qi ^ in order to have epimorphisms M Qi ^ and 
a monomorphism — *■ M — *• (BQi- li I ^ J , \et x G J, x ^ I with ax G I =^ ax = 0, x ^ 0,ya £ a. 
It follows that a is a zerodivisor and thus a G Up, . Also, x ^ li for a certain i otherwise x G fl/j = /. 
But is (co)primary with ax G I C. Ii,x ^ (or ax = 0, .x 7^ in Qi) => a G pi a C pi. In 
particular, if / : m ^ / for m G max{A), then tn = pi for a certain i as before and m G ass{M). 

Q.E.D. 

It is at this precise point that we have to use specific properties of the ring D that will now 
be used in place of A. Indeed, it is well known [4,15,22] that the ring D* = homk{D, k) of formal 
power series is an injcctivc module as there is a canonical isomorphism R = M* = homk{M, k) 2i 
homD^M, D*). From now on and unless specified, we assume that the partial localization has been 
realized. Using therefore k{x) ® M over k{x) in place of M over k, it is thus equivalent to assume 
that M is n-pure, that is dimk{M) < 00 and in this case we have of course M** ~ M and any 
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associated prime ideal is maximum, that is ass{M) C max{D). However, the reader must always 
keep in mind that the original module was pure and thus contained in its localization, that is to 
say no " simplification" is possible in the language of control theory, where we recall once more that 
the transfer matrix approach is superseded by the localization procedure !. 

As wc have soc(M) = (BsoCmiM) where the summation is now done on ass(M) only, in order 
to dualize the short exact sequence soc{M) M ^ M/soc{M) — > 0, we need first dualize 
the various short exact sequences — >■ socm{M) M ^ N ^ 0. However, if 5 is a simple 
module, we have m = ann{S) C ann{S*) C ann{S**) = ann{S). Accordingly, the dual of a simple 
module isomorphic to D/m is an isotypical simple module, because else it would have a proper 
factor module, the dual of which would be a proper submodule of D/m. We get therefore, again 
because of the injectivity of D* , the short exact sequences N* M* — !■ topm{M*) and 
the desired dual sequence is finally obtained by introducing the intersection rad{M*) = CiN* for 
the various m e ass{M). A key result is provided by the following theorem which is not evident 
at all, even on very elementary examples, and provides a link between the socle of a module and 
the top of the corresponding system. 

THEOREM 5.2.2: A^* ~ mM* and the previous short exact sequence relative to m is isomorphic 
to the short exact sequence mM* — » M* D/m ® M* 0. 

Proof. As the second result is just obtained by tensoring with M* the short exact sequence 
O^m^D^D/m^O, it just remains to prove the first one. 

For this, let us set m = (ai,...,at) = X^i'a and use the notations of the preceding lemma. If 
we introduce I : a = J{a) C F and introduce the corresponding short exact sequence 
J{a) ^ F ^ ^{i) ~^ for each generator a G m, we have aJ{a) C / by definition and the 
multiplication by a thus induces a monomorphism — > N{a) M. By duality, we have the 
epimorphism M* A N{a)* and obtain therefore N{a)* = aM*. Finally, if / : m = J C F, 
we have of course J = nJ(a) where the intersection is taken on the various generators of m 
and an exact sequence ^ J ^ F ^ (BN{a) inducing a monomorphism ^ A^ ^ (BN{a) 
because of the defining short exact sequence O^J^F^N^Q. Moreover, the inclusion 
J C J(a),Va G m induces an epimorphism N — > N{a) 0. Accordingly, we are exactly in the 
position to use Proposition 4.7 in order to get by duality the inclusion N{a)* N* and 
therefore A^* = '^N{a)* = '^aM* = mM* C M*. Finally, any nonzero element in the module 
on the right in the sequence of the theorem is killed by m and admits a representative in R which 
is not in mR. By definition of ham, it is the restriction of a section of to a simple submodule of 
M. 

Q.E.D. 

COROLLARY 5.2.3: We have the short exact sequence CimR ^ R top{R) — > coming 
from the Chinese remainder theorem by tensoring £>/ n mj ~ (BD/mi with R. 

COROLLARY 5.2.4: R ®Rm projects onto top{R) = ©topm(i?) = (Btop{R)m. 
Proof. First of all, using the exactness of the localizing functor, we have top{R) = (BD/m (g) R 
with D/m (8) i? ~ R/mR =^ Rm/mRm — {R/mR)m — R/mR ~ topm{R). Indeed, m G max{D) ^ 
(tn, s) = D, Vs e £> - {m} and 3t & D — {m}, a e m with st + a — 1. Accordingly, V/ G R, then 
■j/ = ■^^7^ ^ tf + ^f and we can therefore take out the denominators when localizing. Finally, as 
m' + m" = D,Vm',m" G ass{M), we have similarly {R/m'R)m'^ = 0. 

Q.E.D. 

This corolary allows one to use Nakayama' lemma in order to look for the generators of R 
because M and thus R are finitely generated over k and thus over D by assumption. The following 
theorem, which is a straight consequence of ([12]. IV, §2, p 109), constitutes the secret of Macaulay 
([13], §82, end p 91) and explains the reason for introducing the (inverse) system. 

THEOREM 5.2.5: The minimmn number of generators of R is equal to the maximum number 
of isotypical components in the various components of the socle of M or of the top of i? = M*. 
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5.3) As the examples in this subsection will clearly show out, the number of generators is related to 
the localized module/system and not to the original module/system as we shall exhibit situations 
needing two generators even though max{D) n ass{M) — 0. Therefore, in this last subsection, we 
shall explain the way followed by Macaulay ([13], §79, p 89) in order to get back informations on 
the original system from results on the localized one. For simplicity the index k of the unknowns 
will not be written down. 

As we already did, setting x = (Xi) •••) Xn-r) and introducing = (/z',/i") with 1 < fi'i < 
n — r,n — r + 1 < j < n, we obtain the locahzed system by substituting = y(y^f^") = x^'^m" 
in the original PD equations. However, even if we start from an involutive system, the correspond- 
ing localized system may be still involutive with full classes but with quite different features, for 
example no longer homogeneous if the original system is homogeneous. In order to manage with a 
solved form, we have the following result found by Macaulay ([13], §78, p 88 (A) and §79, p 89 (B)). 

PROPOSITION 5.3.1: The localized system is k{x) ® R- 

Proof. As M is finitely presented, there is the abstract isomorphism ([22], Th 3.84, p 107): 

homk^^){k{x) M,k(x)) ^ Hx) ^ homk{M,k) = kix) (E> R 

where we recall that M can be identified with its image in k{x) ® M as it is pure. However, 
in actual practice, it is not evident at all to discover that a single determinant in k{x) allows to 
provide modular equations with coefficients polynomials in k[x]- In fact, all principal jets (pri) 
of order > q and class > n — r + 1 can be expressed from parametric jets {par) of the original 
system and can therefore be expressed by means of finite linear combinations of the jets of the 
localized system with order > g — 1 with coefiicients in k[x\ {care). However, these latter jets can 
themselves be linearly dependent through a finite number of equations of order < q—1 . Solving 
these equations with respect to principal jets of order < q — 1 may therefore bring a determinant 
c(x) G k\x\- Accordingly, any modular equation of the localized system can be written in the form 
E = c{x)a'P"''^°- + ^h{x)aP" = with 6, c £ k[x\- More generally, we get a finite number of modular 
equations of the form E = c^^^{x)a^' = c^„x\'a^" = called r-dim,ensional modular equations by 
Macaulay, with an inequality | A' | — | /i" |< 5 for a certain relative integer 5 and equality for ho- 
mogeneous systems. As no ^^simplification!'' may exist, that is tr-i{M) = for the original module, 
one just needs to set a'*" = Xm'"'^ '^ ^ order to get the so-called n- dimensional modular equa- 
tions Ea' = X]A'+^i'=a'''M""^'^ ''^"^ ~ ^" following theorem on the way to generate the modular 
equations is the key result obtained by Macaulay ([13], §82, end p 91) but not understood up to now. 

EXAMPLE 5.3.2: Purity is essential in the process of localization/ delocalization. In order to 
prove it, let us consider the very simple Example 2.3 of codimension 1 but not pure. Localizing, 
we get from the second PD equation XiV'i = leading to 7/2 = that is a new PD equation. As an 
exercise, we let the reader compare with the pure situation of Example 2.2. 

THEOREM 5.3.3: There is a finite number of r-dimensional modular equations, a smaller num- 
ber of r-dimensional modular equations of which all the others are derivates and an equal or still 
smaller number of n-dimensional modular equations of which all the others of a certain order are 
derivates. 

Proof In we have | |=| /x' | -f | /i" | = | a' | - | A' | -|- | /i" |>| a' \ -6. As for the derivates, 
we have dyE = c^Ixa'+m'^^"'''*""^"^ = with a''-^' = if /Zi = or a(''i'->'''-i'-''''') ii m > I 
([14], §60, p 69). 

We have therefore {dyE)/}, = dY'{Ef3') = dyE^, = Ea'-f/^'=/3'C^''«^''' = with [7" |< r 
since there are only a finite number of linearly independent derivates of the r-dimensional equations 
as the localized system is a finite dimensional differential vector space over k{x)- 
More precisely, starting now from Ea' = c^,^a^°' '^"^ = with A' < a', the 7"-derivates are of 

the form Ef}, = cKa^/^'-^'^i^" = with A' < /3',7" < /x" and where | /3' |, | 7' |< t are fixed 
positive integers (including zero). 

Let us finally consider all the modular equations of order q that can be obtained as derivates, 
that is all the E^, with g >| /3' | - | A' | | /x" | - | 7^' |>| /3' | -t - S that is to say 
P'i <| |< q + S + T,yi = 1, n — r. Accordingly, every E/}i of order g is a derivate of E^' for a 
fixed a' if a'i> q + S + T,\/i = 1, n — r. 
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Q.E.D. 



EXAMPLE 5.3.4: Coming back to Example 5.1.2, we have par = {y, yi, y2, J/ii} ^ M ~ 
ky + kyi + ky2 + kyu and thus / = (1, 0, Q,{)) ^ Ei = = f = (0, 1, 0, {)) -> E2 = = 
0, / = (0, 0, 1, 0) ^ = a2 = 0, / = (0, 0, 0, l)^ Ei = a" = 0. We have soc{M) = Dy2 © Dyn 
with two isotypical components both killed by m = {di,d2) and top{R) = {i?3,i?4} provides two 
generators for R as we have indeed 1^2 i^s = Ei,diE4 = E2, that is to say mR is generated by 
{El, E2} in agrement with Nakayama's lemma. 

EXAMPLE 5.3.5: Coming back to Example 5.1.3, we have par ~ {y, yi, 2/2, 2/3, J/ii, 2/12, 2/13, J/iii}- 
We have soc{M) = Dym killed by m = (^1,^2) and top{R) = {E} with E = a^^^ + a^^^ = 
provides a unique generator of R as we have diE = a^^ + a^^ = 0, d2E = a^^ = 0, rfaE' = a^^ = 
0, ...,diiiE = a° = as a way to generate mi?. It is remarkable that 8 = 2^ = 2" is a general 
combinatorial result proved by Macaulay ([13], §68, p 79). A similar simpler situation is met with 
n = 2 and y22 = 0, yi2 — yii =0 (exercise). 

EXAMPLE 5.3.6: Coming back to Example 3.14 which needs a partial localization with k{x) = 
k{xi,X2,X3), the localized system is y| = 0, y| = 0, y| = 0,X3y^ + X2y^ + XW^ = 0- Clearly 
k{x) ® M ~ k{x)y^ + k{x)y^ — k{x)[di\y^ © k{x)[d4\y^ is a semi-simple module with two iso- 
typical components both killed by m = (^4). Accordingly, (1,0) — > xs-^i = Xi'A ~ = 
0, (0, 1) X3E2 = X30-2 ^ X2a3 = provides the two generators of the locaUzed system, even 
though max{D) n ass{M) = in this case. We notice that the determinant c(x) = X3 is unavoid- 
able. Delocalizing, we get X1X3 a} — a| = 0, X2X3 ^ — 03 = and so on, in agrement with 
the general theory for the original system. 

EXAMPLE 5.3.7: n = 3,m = l,q = 2,k = Q, the module is 2-pure and defined by the homoge- 
neous involutive system y^s = 0, j/23 — yi3 = 0, 1/22 — yi2 = 0. With k{x) = k{xi ) , the corresponding 
localized system y33 = 0,y23 - Xi^a = 0,y22 - Xil/2 = is again involutive with par = {y, 2/2, 2/3}- 
We obain therefore (1,0,0) Ei = a° = 0,(0,1,0) ^ E2 = a"^ + xia^^ + (xi)^a222 ^ ^ 
0,(0,0,1) ^ E3 = a^ + xia23 -|- (xi)^a^^^ + ... = 0. We notice that mi = (^3,^2 - xi) kills 
ys while m2 = (^3,^2) kills y2 — Xi?/, each max:imum ideal in fc(xi)[c?2, <^3] leading to a unique 
isotypical component. Denoting simply by M the localized module and by R the corresponding 
system, we have the short exact sequence — > socm2{M) ^ M ^ N2 ^ (care to the notation) 
and the dualizing short exact sequence R ^ topxa^^R) — * 0. Here, N2 is obtained 

by adding y2 — Xi2/ = to the equations of the localized system and we obtain the subsystem 
m = foiEi + X1E2) + fsEs CR = foEi + /2^2 + /3^3. We check the relations: 

d2-Bi = 0, d^Ei = 0, d2E2 = Ei + xiE2,d3E2 = 0, ^2-^3 = XiE3,d3E3 = Ei + X1E2 

transforming R into a 3-dimensional differential vector space over k{x) and obtain m2i? = (/2 + 
/3)(i^i +Xi^^2) +Xi/3i^3 = f'oiEi+xiE2) + .f'3E3, that is wc check directly = and could 
check similarly = miR, a result highly not evident at first sight. According to the general 
theory, there should be one generator only and we may choose E = E2 + E3 in order to generate 
R as we have indeed the three linearly independent relations: 

E = E2 -\- E3, d2E — xiE = El, dsE = Ei -\- X1E2 

allowing to determine Ei,E2, E3 from the derivates of E. The system being homogeneous, we have 
g + (5 + T = 2 - 1 + 1 = 2. As wc have £; = l(ai + a3) + ... + {xif{a^^^^ + a^^^ + 0^22 + + ^123 + 
a223) + ... = 0, it is easy to check that the single modular equation En = a^^2 _|_ ^122 _|_ ^222 _|_ 
Qjii3 ^ ^123 ^ ^223 ^ generates o^^ = 0, o^^ -|- a'^^ = 0, + ^22 ^ ^1 ^ ^2 ^ q ^^3 ^ q ^^0 ^ q 
successively. Hence all the modular equations up to order 2 are generated by a single modular 
equation at order 3, a result highly not evident at first sight. 

EXAMPLE 5.3.8: Looking back to Example 2.3, the primary decomposition brings the two 
subsystems R' = {y2 = ^ a° = 0,a^ = 0, a" = 0,...} and R" = {y22 = 0,yi2 = 0,yii = ^ 
a° = 0, = 0, a2 = 0} with R = R' + R" . One clearly needs two generators in order to generate 
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any Rq, say {o^ = 0, o^^^ = 0} for q = 3. 

EXERCISE 5.3.9: Treat similarly the inhomogeneous involutive system j/33 = 0,t/23 = 0,t/22 = 
0, 2/13 — 2/2 = and prove that there is again a unique generating modular equation. 

6) CONCLUSION: 

We summarize the way leading to revisit the "inverse system" of Macaulay by using modern 
techniques of algebraic analysis, namely differential geometric arguments for studying the system 
instead of the module. 

The main purpose is to find generators for the differential system dual to the differential module. 
For this, the only way known in the literature is to control the generators of the system from the 
generators of its top by using Nakayama's lemma. Again by duality, this amounts to count the 
isotypical components of the socle of the module. Meanwhile, the key idea is to decompose the 
system into subsystems instead of using a primary decomposition of the module in order to deal 
with pure modules, a concept generalizing the unmixedness assumption of Macaulay. However, 
the original system is not in general finitely generated and it is therefore essential to use a partial 
localization in order to deal with a finite dimensional localized system. 

The present approach avoids the abstract systematic use of the injective hull by Oberst and opens 
a new way towards effective computer algebra packages for studying identifiability in engineering 
sciences. 

It is thus remarkable that Macaulay had the intuition of these techniques as early as in 1916 and 
we express our deep gratitude to his work. 
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